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Abstract. For High frequency BL Lac objects (HBLs) like H1426+428, a significant fraction of their TeV gamma- 
rays emitted are likely to be absorbed in interactions with the diffuse IR background, yielding e pairs. The 
resulting e pairs generate one hitherto undiscovered GeV emission by inverse Compton scattering with the 
cosmic microwave background photons (CMBPs). We study such emission by taking the 1998-2000 CAT data, 
the reanalyzed 1999 & 2000 HEGRA data and the corresponding intrinsic spectra proposed by Aharonian et al. 
(2003a). We numerically calculate the scattered photon spectra for different intergalactic magnetic field (IGMF) 
strengths. If the IGMF is about 10 _18 G or weaker, there comes very strong GeV emission, whose flux is far above 
the detection sensitivity of the upcoming satellite GLAST! Considered its relatively high redshift (z = 0.129), 
the detected GeV emission in turn provides us a valuable chance to calibrate the poor known spectral energy 
distribution of the intergalactic infrared background, or provides us some reliable constraints on the poorly known 
IGMF strength. 
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1. Introduction 

Several High frequency BL Lac objects (HBLs) such 
as Mkn501, Mkn421, PKS 2155-304, 1ES 2344+514, 
; H1426+428 and 1ES1959+650 are of great interest, 
, since they emit TeV photons (Catanese & Weekes 
1999;Aharonian et al. 1999, 2002a; Horns 2002). As early 
as in 1960s, it was pointed out that the observed high 
energy gamma photon spectra from TeV sources might 
be significantly modified — as the original high energy 
gamma-rays satisfying E^Eb > 2(m c c 2 ) 2 (where £J 7 is 
the energy of the seed high energy gamma-ray, Eb is the 
infrared— UV background photon) travel through the in- 
tergalactic infrared— UV background, a significant fraction 
of them will be absorbed, leading to e^- pairs (Nikishov 
1962; Gould & Schreder 1965; see also Stecker, de Jager 
& Salamon 1992). If the optical depth to TeV photons 
is evaluated, the intrinsic spectrum can be derived. Such 
calculations were firstly made for Mkn501 (during the 
1997 flaring activity), and then for H1426+428, indicat- 
ing that their intrinsic hight-energy spectra have broad, 
flat peaks which are much higher than the observed one 
in the ~ 5 — 10 TeV range (e.g. Konopelko et al. 1999; 
Aharonian, Timokhin & Plyashcshnikov 2002; de Jager 
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& Stecker 2002; Aharonian et al. 2003a). Very recently, 
Dai et al. (2002) have suggested that inverse Compton 
(IC) scattering of the resulting pairs against the cosmic 
microwave background photons (CMBPs) may produce a 
new GeV emission component in TeV blazars . As an 
illustration, they have calculated the well-studied blazar 
Mkn501, and obtained strong GeV emission, as long as 
the intergalactic magnetic field (IGMF) is weak enough. 

In this paper, we turn to investigate another impor- 
tant extragalactic object, H1426+428, which has been de- 
tected by Remillard et al. (1989) firstly. With the rela- 
tively high redshift z = 0.129, the optical depth exceeds 

For Gamma-ray bursts, such Compton scattering leads to 
an observable, delayed MeV-GeV emission component. In fact, 
there are more than six events have been detected, the well 
known one is the GRB 940217, whose delayed MeV-GeV emis- 
sion lasts for ~ 5000s (Hurley et al. 1994). If the model pro- 
posed by Plaga (1995) and Cheng & Cheng (1996) (see also 
Dai & Lu 2002; Wang et al. 2003) is correct, the intergalactic 
magnetic field is « 10" 20 G (For GRB 940217, the redshift is 
unavailable. But with the empirical relation established based 
on several GRBs with known redshifts, a rough estimate sug- 
gests its redshift is ~ 0.68 (Fenimore & Rumirez-Ruiz 2000). 
On the other hand, the typical energy of the observed delayed 
photons is 50MeV, which implies the typical energy for the 
high energy seed photons E-, ~ 0.3 TeV. Substituting these 
two estimates into equation (3), we have Big ~ 10 _20 G). 
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unity for energies of ~ 300GeV (e.g. de Jager & Stecker 
2002). Any detection of a signal at TeV energies translates 
directly into a high luminosity of the source. The detec- 
tion and spectral measurement of H1426+428 by the CAT, 
and VERITAS groups (Djanati-Atai et al. 2002; Petry et 
al. 2002) indicate a steep spectrum between 250GeV and 
ITeV, whereas at higher energies, detections have been 
reported by the HEGRA group based upon observations 
carried out in 1999 and 2000 (Aharonian et al. 2002b). 
Here, following Dai et al. (2002), we study such hitherto 
undiscovered GeV emission by taking the 1998-2000 CAT 
data (Djannati-Ata'i et al. 2002), the reanalyzed 1999 & 
2000 HEGRA data as well as the derived intrinsic spectra 
proposed by Aharonian et al. (2003a). 

The predicted GeV emission is sensitive to IGMF 
(Dai et al. 2002). If the IGMF is stronger than 10~ 12 G, 
the formation of a very extended electron/positron Halo 
(i? ~ lOMpc for the primary photons with a typical en- 
ergy ~ ITeV) is unavoidable (Aharonian, Coppi & V61k 
1994), and its emission is nearly isotropic. If this is the 
case, the GeV emission predicted below will be signifi- 
cantly lowered. However, the strength of IGMF has not 
been determined so far. Faraday rotation measures imply 
an upper limit of ~ 10~ 9 G for a field with 1 Mpc correla- 
tion length (sec Kronberg 1994 for a review). In dynamo 
theories, to interpret the observed ^G magnetic fields in 
galaxies and X-ray clusters, the seed fields required may 
be as low as 10~ 20 G (see Kulsrud 1999 for a review). 
Theoretical calculations of primordial magnetic fields gen- 
erated during the cosmological QCD or clectroweak phase 
transition show that these fields could be of order 10~ 20 
G or even as low as 10 -29 G (Sigl, Olinto & Jedamzik 
1997). The model proposed by Aharonian, Timokhin & 
Plyasheshnikov (2002) to explain the observed TeV emis- 
sion of Mkn501 suggests that B\q < 10~ 18 G. Modelling 
the long delayed MeV-GeV emission of GRB 940217 sug- 
gests that Big ~ 10 _20 G (see footnote 1). It is evident 
that all of these claims are far from conclusive. Therefore, 
observing a hitherto undiscovered GeV emission compo- 
nent from flares of TeV blazars such as H1426+428, one 
may be able to obtain important information or con- 
straints on the poorly known IGMFs. 

This paper is structured as follows: In section 2, wc 
describe our physical model briefly. In section 3 we provide 
our numerical results. In section 4 we give our conclusions. 

2. Physical model 

As mentioned above, as the 7— rays with the energy high 
up to ~ ITeV travel toward the observer, a significant 
fraction of them is likely to be absorbed, yielding e ± . 
The pair production optical depth depends strongly 
on the 7-ray energy 2 (E' = (1 + z)E y _ b s ) and the red- 
shift (z) (throughout this letter, the subscripts "obs" and 

2 In the previous works (e.g. Dai et al. 2002), the red-shift 
correction of j c as well as the four timescales involved (see 
below) has been ignored. However, in the present case, for the 



"loc" denote the parameters measured in the observer 
and source frames respectively) . If a primary photon with 
energy E' has been absorbed, the resulting pairs 
have Lorentz factors 7 e ~ ^/(2m e c 2 ) w 10 6 £^/1 TeV), 
where m e is the electron mass. Such ultra-relativistic 
pairs will subsequently Compton scatter on the ambient 
CMBPs, and boost them to an average value ~ 7 2 e ~ 
0.63(1 +z)(E! y /l TeV) 2 GeV, where e = 2.7fcT is the mean 
energy of the CMB photons with T ~ 2.73(1 + z)K and k 
is the Boltzmann constant. 

Since the emission of blazars is jet-like rather than 
isotropic, so do the scattered CMBPs. Thus it is neces- 
sary to investigate whether significant part of them has 
been scattered out of the line of sight (the "escaping" ef- 
fect), or not. For the IC- scattered CMBPs, which move 
nearly along the path of the ultra-relativistic electrons and 
the deflection angle is <~ I/70, much less than the typical 
open angles of blazars. On the other hand, the recoiling 
of the seed electrons is small, too. It is easy to show that 
for one scattering, the deflection angle is < e/m c c 2 . For a 
seed electron scattering N times, its final deflection angle 
can be estimated by 9d ~ VNe/m e c 2 . Typically, for the 
IC scattering considered in this paper, N rts 1000, which 
implies #d ~ 10~ 7 . Therefore, the "escaping" effect can 
be ignored. 

2.1. The GeV emission duration 

As shown in Dai et al. (2002), there are four timescales 
involved in the emission process: The first is i V ar.obs 7 the 
observed variability time of the source emission, which has 
been mentioned earlier. 

The second is the well-known angular spreading time 

Ai A , obs « (l+z)%!£ = 960(l+z)(-^)- 2 (-^)- 1 s,(l 
27^0 10° 0.1cm J 

where R pail = (0.26ctt^ir) 1 « 5.8 x 10 25 (^a^ r )- 1 cm 
is the typical pair-production distance, <tt is the Thomson 
cross section, and niR, ~ 0.1cm -3 is the intergalactic in- 
frared photon number density (Dai & Lu 2002; Dai et al. 
2002; Wang et al. 2003). 

The third is the IC cooling timescale. In the absence of 
any IGMF, the IC cooling timescale in the observer frame 
would be (Dai et al. 2002) 

At IC , obs ~ (l+z)t ICJoc /(2 7 2 ) = 38(l+z)- 3 ( 7e /10 6 )- 3 s.(2) 

tiC,ioc 5 the IC cooling time scale measued in the source 
frame is estimated by iic.ioc = 3m e c/(47 e <7T'u cm b) = 7.7 x 
10 13 (1 + z)- 4 (7 e /10 6 )" 1 s, where u cmb = aT 4 is the CMB 
energy density, and a is the radiation constant. 

The fourth is the magnetic deflection time, which ap- 
pears in the presence of IGMFs. The additional emission 
time due to magnetic deflection is (Plaga 1995; Dai et al. 
2002) 

A^B.obs (l/2)(l + *)ti C 

relatively high red-shift of H1426+428, z = 0.129, the red-shift 
correction must be taken into account. 
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for #b, the deflection angle ~ Aic/i?L = 1.3 x 
10- 5 ( 7e /10 6 )- 2 (Bi G /10- 20 G)(l + z)" 4 < 1, where Ai C ~ 
ctic.ioc, Rl — 7e"ieC 2 /eSiG is the Larmor radius of the 
electrons, Big is the strength of the IGMFs. 

Consequently, the duration estimate of the IC emis- 
sion from electron/positron pairs scattering off the 
CMB is Ai(7 e ) = max(At IC , bs, At A ,obs, Ai B ,obs, ivar.obs)- 
Generally, the typical duration of GeV emission is de- 
termined by the source activity time and the energy de- 
pendent magnetic deflection time. For convenience, we in- 
troduce a critical Lorentz factor 7 CjC which is defined by 

AtB.obs ^var.obs 



7c, c 

Toe 



0.59(%^)" 
lday 



1/5 



Big 
10- 20 G 



) 2 ' b {l + z) 



-11/5 



(4) 



For electrons with Lorentz factor 7c < 7 C . C , the magnetic 
deflection angle is such large that most of the scattered 
CMBPs can not reach the observer in a finite time t V ar,obs 7 
which contribute little to the GeV emission considered 
here. Therefore, we take 7e, c as the lower limit for the 
integration performed below. 

2.2. Basic formulae for calculation 

For the ultra-relativistic seed electrons with the distribu- 
tion dN e /d'y e , the time-averaged scattered photon spec- 
trum is given by (Blumenthal & Gould 1970; also see Dai 
& Lu 2002; Dai et al. 2002) 



dN: 



sc 



dEy 



1 



V dj e J \dtdEy J 



dtdj c 



(5) 



where is the luminosity distance to the source, E~y> is 
the externally scattered photon energy, t is the time mea- 
sured in the local rest frame, and (Blumenthal & Gould 
1970) 



dN 7e,e nrlcn(e}de E Y 
dtdEy ~ ^^H 2 ^' 111 ^) 



k 4 7e V 

+E Y +4 7 2 e -B 2 /2 7c 2 6], 



is the spectrum of photons scattered (using the Tompson 
cross-section formula) by an electron with Lorentz factor 
of 7e from a segment of the CMBP gas of differential num- 
ber density (Blumenthal & Gould 1970) 



n(e) = [ 7 r 2 (n C ) 3 ]- 1 [e 2 /(e e 



!)]• 



(7) 



The integration for variable 7 ranges from 7cjC to 



K PiCUt /2m c 2 



where E' 



up t is the observed upper cut- 
off energy timed by (1+z). Please note that this is only 
valid at the beginning since j c decrease with t, which 
can be written into 7 e ( t ) = 7e ,o/(l + &7e,oi)i where 7c ,o 
is the original Lorentz factor of the seed e* 1 pairs, b = 
4a T U cmh /3m c c. 

Dai & Lu (2002) simply replace J dt by iic,ioc ; and 
assume the seed electrons have a distribution oc 7<T ri . 



With these simplifications, they have shown that in 
the Thomson limit, the scattered photon spectrum oc 



E ^T 1+ 2)/2 for A £ B;Qbs ( 7e ) <<; tyai ohs and K E JPi 



-3)/2 



for AtB,obs(7o) ^ ivar,obs- Instead of taking such a sim- 
ple estimate, below we make a detailed numerical cal- 
culation: the integration for variable t ranges from to 
7.7 x 10 13 (1 + z)- 4 ( 7 o !C /10 6 )^ 1 s . 



3. Numerical result 

In our calculation, the original seed electron energy spec- 
trum is derived as follows: we simply apply a polynomial 
fit to the local highly absorbed photon spectrum: Setting 
x = log(^/1.0TeV), we have 



dN 7 
dE^ 



10 



0.79046-3. 31624a; + 1.15145z -0.79715x 



(8) 



which is in units of AirD^lO 



12 ph cm 2 s 

local intrinsic photon spectra dN* int /dEy 



TeV" 1 . The 
(also in units 

ph cm -2 s _1 TeV -1 ) can be approximated 
as follows (Aharonian et al. 2003a): 



of 4ttL>£10- 12 



dN^ 



EL 





■■ 8xl *W>~ 


dE^ ~ 5 


^int „ 1 


[ i. 6x io 2 (o | fev 


dE' y ~ | 




dN iMt „ 1 


f 6 .5 x lO^oi* 


dE^ ~ 1 


I 2 ^- 4 (l.60Tov)' 



■3/2 



(9) 



E; < 0.85TcV, 
- 2 , E; > 0.85TeV. 



E; < 1.6TcV, 
E' > 1.6TV. 



where the superscript i = a,b,c represent the correc- 
tion models proposed by Primack et al. (2001), Aharonian 
et al. (2003a) and Malkan & Stecker (2001) respectively. 

l(7 e /10 6 ) TeV, we have (including positrons) 



Since 



Ml 

dj e 



2 x 10~ 6 ( 



<^,int 

dEL 



dN-, 
~dE~i 



(12) 



(6) Now we can calculate the spectra of the scattered CMBPs 



with equations (5)-(12). Please note that the e, Ey in 
equations (6), (7) are all measured in the local frame 
rather than by the observer. 

The numerical result has been shown in Figure 1. 
As long as the IGMF strength is weak enough, there is 
the very strong GeV emission (although it is intrinsic 
spectrum-dependent), whose flux is far above the thick 
dashed line, the GLAST sensitivity computed for an expo- 
sure time of 2 days. Just as shown in figure 1, the inferred 
intrinsic spectra is most sensitive to the characteristics 
of the extragalactic background light. As a consequence, 
different correction models lead to much different intrin- 
sic spectra, so does the accompanying GeV emission, and 
thus the detection of such emission may in turn provide us 
a valuable chance to test those IR background absorption 
models. 
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Fig. 1. The time-averaged high energy Gamma-ray spectra of 
H1426+428 for a flare with t V ar,ob s = 2 days. The filled cir- 
cles are the reanalyzed f999/2000 spectral data measured by 
HEGRA ( Aharonian et al. 2003a) , and the filled rectangles are 
the f998-2000 spectral data measured by the CAT (Djannati- 
Atai' et al. 2002). In a, b and c: the thick solid lines represent 
the derived intrinsic spectra based on the correction model 
proposed by Primack et al. (2001), Aharonian et al. (2003a) 
and Malkan & Stecker (2001) (see also figure 4 of Aharonian 
et al 2003a) respectively; the thick dashed lines represent the 
GLAST sensitivity computed for an exposure time of 2 days; 
the dotted lines are the secondary photon spectra calculated 
in the text from the resulting e pairs interacting with the 
CMBPs, the number 1 — 3 correspond to IGMF strength of 
10~ 20 G, 10~ 19 G and 10~ 18 G, respectively. 

In our calculations, as an illustration, we take 
-^up cut = 10-2 (1 + z)TcV (Above which there is no obser- 
vation available), which is conservative. Even so, with the 
three correction models taken in Aharonian et al. (2003a), 
the predicted time-averaged scattered photon energy flux 
are far above (case b, c), or at least comparable with (case 
a) that of the well studied Blazar Mkn501 (see Dai et al. 
2002 for detail). One problem arises, just as shown in fig- 
ure 1 c, the predicted new hard 7— ray emission at the en- 
ergy ~ 0.24 TeV is about ten times that observed. In prin- 
ciple, such emission may be absorbed by the background 
IR/UV photons again, but the resulting flux is still high 
above the observed (for E 1 ~ 0.24TeV, - 1). Two 
possibilities exist: one is that the derived intrinsic spec- 
trum based on the correction model proposed by Malken 
& Stecker (2001) is much harder than what it really is. 
Another is that the IGMF might be strong enough, the 
observable GeV- TeV emission flux has been suppressed 
significantly. 

For Blazar H1426+428, due to the poorly known X-ray 
band spectrum, we cannot provide a definite SSC spec- 
trum as Dai et al. (2002) have done for Mkn501. Even 



so, as shown in Figure 1, at a first glance, the SSC com- 
ponent is likely far below than the externally scattered 
photon flux. 

4. Conclusions 

The HBLs HI426+428 is distinguished by its relatively 
high rcdshift as well as its strong TeV energy emission. 
In this Letter, with the recently observed data as well as 
the derived intrinsic spectrum, we have studied the pos- 
sible accompanying GeV emission, which is due to the IC 
scattering of CMBPs by the pairs produced in inter- 
actions of high energy photons with the cosmic infrared- 
UV background photons. If the IGMF is weak enough, 
there is very strong hitherto undiscovered GeV- TeV emis- 
sion, whose flux is far above or at least comparable with 
that of the well-studied HBLs Mkn501 (Dai et al. 2002). 
With the upcoming satellite GLAST, planned for launch 
in 2005, it is easy for us to detect such emission directly. 
We also noted that the shape of the inferred source spectra 
is most sensitive to the characteristics of the extragalac- 
tic background light. A different correction model leads to 
a different intrinsic spectrum, so does the accompanying 
GeV emission (see figure 4 for detail). Thus the detection 
of such emission may in turn provide us a valuable chance 
to test those IR background absorption models. 

In Aharonian et al. (2003b), the detection of TeV 
7— rays from the BL Lac 1ES1959+650 has been reported. 
1ES1959+650 is located at a redshift of z = 0.047, provid- 
ing an intermediate distance between the nearby blazars 
Mkn42I and Mkn504, and the much more distant object 
H1426+428. Without doubt, for such a source, the ob- 
served TeV emission has been absorbed significantly by 
the infrared- UV background photons, and thus a strong 
GeV emission is expected, as long as the IGMF is weak 
enough. 
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